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Abstract 

a 

We extend Goldie's (1991) Implicit Renewal Theorem to enable the analysis 
<^ of recursions on weighted branching trees. We illustrate the developed method 

I by deriving the power tail asymptotics of the distributions of the solutions R 

to 

JV / JV \ 

R = Y, CiRi + Q, R = ( V CiRi v Q ' 

i=l \i=l / 

and similar recursions, where (Q, N, Ci, C2, . . . ) is a nonnegative random vector 
with N £ {0, 1, 2, 3, . . . } U {00}, and {i?j} ig N are iid copies of R, independent 
of (Q, N, Ci, C2, • • • ); here V denotes the maximum operator. 
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1. Introduction 

(N 

C*~) This paper is motivated by the study of the nonhomogeneous linear recursion 

R = ^2CiRi + Q, (1) 

where (Q, JV, C\, C2, ■ ■ ■ ) is a nonnegative random vector with jV £ NU {00}, N = 
{0,1,2,3,...}, P(Q > 0) > 0, and {i?i}i S N is a sequence of iid random variables, 
independent of (Q, JV, C\, C2, . . . ), having the same distribution as R. This recursion 
appeared recently in the stochastic analysis of Google's PageRank algorithm, see [2"7l 
HU] and the references therein for the latest work in the area. These types of weighted 
recursions, also studied in the literature on weighted branching processes [25] and 
branching random walks [8] , are found in the probabilistic analysis of other algorithms 
as well e.g., Quicksort algorithm [13]. 

In order to study the preceding recursion in its full generality we extend the implicit 
renewal theory of Goldie [TJ to cover recursions on trees. The extension of Goldie's 
theorem is presented in Theorem |3.1| of Section [3] One of the observations that allows 
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this extension is that an appropriately constructed measure on a weighted branching 
tree is a renewal measure, see Lemma 3.1 and equation In the remainder of the 



paper we apply the newly developed framework to analyze a number of linear and 
non- linear stochastic recursions on trees, starting with ([I]). Note that the majority of 
the work in the rest of the paper goes into the application of the main theorem to 
specific problems. 

In this regard, in Section |4j we first construct an explicit solution (17) to ([!]) on 



a weighted branching tree and then provide sufficient conditions for the finiteness 



of moments and the uniqueness of this solution in Lemmas 4.4 and 4.5 respectively. 
Furthermore, it is worth noting that our moment estimates are explicit, see Lemma |4.3[ 
which may be of independent interest. Then, the main result, which characterizes the 
power-tail behavior of R is presented in Theorem |4.1| In addition, for integer power 
exponent (a £ {1, 2, 3, . . . }) the asymptotic tail behavior can be explicitly computed as 
stated in Corollary |4.1| Furthermore, for non integer a, Lemma |4. 1| yields an explicit 
bound on the tail behavior of R. Related work in the literature of weighted branching 
processes (WBPs) for the case when N = oo and Q, {Ci} are nonnegative deterministic 
constants can be found in I25j (see Theorem 5), and more recently, for real valued 
constants, in [5]. However, these deterministic assumptions fall outside of the scope of 
this paper; for more details see the remarks after Theorem |4.1| in Section |4~2) 

Next, we show how our technique can be applied to study the tail asymptotics of 
the solution to the critical, E YliLi Q =1, homogeneous linear equation 

N 

R = Y,CiRi, (2) 

i=l 

where (N, C\, C%, . . . ) is a nonnegative random vector with N € NU{oo} and is 
a sequence of iid random variables independent of (N, C%, C2, . . . ) having the same dis- 
tribution as R. This type of recursion has been studied to a great extent under a variety 
of names, including branching random walks and multiplicative cascades. Our work is 
more closely related to the results of [53] and [T7], where the conditions for power-tail 
asymptotics of the distribution of R with power exponent a > 1 were derived. In 
Theorem |4.2| of Section |4.2| we provide an alternative derivation of Theorem 2.2 in 
[23] and Proposition 7 in |17j . Furthermore, we note that our method yields a more 
explicit characterization of the power-tail proportionality constant, see Corollary |4.2[ 
For the full description of the set of solutions to |2| see the very recent work in [3J. 
For additional references on weighted branching processes and multiplicative cascades 
see [2 [23l [22l [2SJ [24] and the references therein. For earlier historical references see 

[201 na mi. 

As an additional illustration of the newly developed framework, in Section [5] we 
study the recursion 

r£ njCiluj V Q, (3) 

where (Q, N, C\, C2, ■ ■ ■ ) is a nonnegative random vector with N E N U {00}, 
P (Q > 0) > and is a sequence of iid random variables independent of 

(Q, N, C\, C2, . . . ) having the same distribution as R. We characterize the tail behavior 
of P(R > x) in Theorem |5.1 [ Similarly to the homogeneous linear case, this recursion 
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was previously studied in |B] under the assumption that Q = 0, N = oo, and the 
{C,} are real valued deterministic constants. The more closely related case of Q = 
and {Ci} > being random was studied earlier in |18) . Furthermore, these max-type 
stochastic recursions appear in a wide variety of applications, ranging from the average 
case analysis of algorithms to statistical physics; see p] for a recent survey. 

We conclude the paper with a brief discussion of other non-linear recursions that 
could be studied using the developed techniques, including the solution to 



First we construct a random tree T. We use the notation to denote the root 
node of 7", and A n , n > 0, to denote the set of all individuals in the nth generation 
of T, Aq = {0}. Let Z n be the number of individuals in the nth generation, that is, 
Z n = \A n \, where | • | denotes the cardinality of a set; in particular, Zq = \. 

Next, let N + = {1, 2, 3, ... } be the set of positive integers and let U — UfcLo(^+) fc be 
the set of all finite sequences i = (ii, %2, ■ ■ ■ , i n ), where by convention = {0} contains 
the null sequence 0. To ease the exposition, for a sequence i = (ii,?2, • • ■ ,ik) € U we 
write i|n = (ii, i%, . . . , i n ), provided k > n, and i|0 = to denote the index truncation 
at level n, n > 0. Also, for i 6 A\ we simply use the notation i = ii, that is, without 
the parenthesis. Similarly, for i = (ii,..-,i n ) we will use = (ii, ■ ■ ■ , i n , j) to 

denote the index concatenation operation, if i = 0, then = j. 

We iteratively construct the tree as follows. Let N be the number of individuals 
born to the root node 0, N$ — N, and let {A^jigc/ be iid copies of N. Define now 

Ax = {i e N+ : 1 < i < N}, 

A n = ■ ■ ■ ,i n ) € U : (H,...,i n _i) G A n _i,l < i n < iV^,.. .,,„_!)}■ (4) 

It follows that the number of individuals Z n = \A n \ in the nth generation, n > 1, 
satisfies the branching recursion 



Now, we construct the weighted branching tree 7q,c as follows. The root node is 
assigned a vector (Qg, iVg, C($^, C(0,2), ■ ■ ■ ) = (Q,N, Ci,C 2 , ■■■) with N e N U {oo} 
and P(Q > 0) > 0; N determines the number of nodes in the first generation of 
T according to Q. Each node in the first generation is then assigned an iid copy 
(Qi, Ni t C(i,i),C(i,2)> ■ • ■ ) of the root vector and the {iVj} are used to define the second 
generation in T according to Q. In general, for n > 2, to each node i € A n -%, we 
assign an iid copy (Qi, N\, Ca,i), C(i,2)) • • • ) of the root vector and construct A n — 
{(i, i n ) : i G A n _i,l < i n < Ni}; the vectors (Qi, Ni, C^i), C( i)2 ), • • • ), i € A n -\ are 
chosen independently of all the previously assigned vectors (Qj, Nj, Crjn, Cn^), • • ■ ), 
j G Ak,0 < k < n — 2. For each node in Tq,c we also define the weight !!(»!,...,»„) via 



\i=l / 

The majority of the proofs are postponed to Section [7] 




2. Model description 



Z n = Ni - 
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the recursion 

n n = c n , n^...^) = C(j l! ... ) j ra )ii(j li ... i j n _ 1 ), n > 2, 

where II = 1 is the weight of the root node. Note that the weight II(j 1) ...^ n ) is equal 
to the product of all the weights C(.) along the branch leading to node (ii, . . . ,i n ), 
as depicted in Figure [TJ In some places, e.g. in the following section, the value of Q 
may be of no importance, and thus we will consider a weighted branching tree defined 
by the smaller vector (TV, Ci, C2, . . . ). This tree can be obtained form Tq.c by simply 
disregarding the values for and is denoted by Tc- 




Figure 1: Weighted branching tree 



Studying the tail behavior of the solutions to recursions and fixed point equations 
embedded in this weighted branching tree is the objective of this paper. 



3. Implicit renewal theorem on trees 

In this section we present an extension of Goldie's Implicit Renewal Theorem [TJ] 
to weighted branching trees. The observation that facilitates this generalization is the 
following lemma which shows that a certain measure on a tree is actually a product 
measure; a similar measure was used in a different context in [§] . Its proof is given in 
Section [TIT] for completeness. Throughout the paper we use the standard convention 
Q logO = for all a > 0. 

Lemma 3.1. Let Tc be the weighted branching tree defined by the nonnegative vector 
(N, Ci, C 2 , . . . ), where iVeNU {00}. For any n€N and i G A n , let V" ; = log IF. For 
a > define the measure 



Hn{dt) = e at E 



ieA n 



n = l,2, 



and let rj(dt) = j-ii(dt). Suppose that there exists j > 1 with P(N > j, Cj > 0) > 
such that the measure P(log Cj G du, Cj > 0, N > j) is nonarithmetic, 
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< E 



< oo and E 



HiLi CI = 1- Then, rj(-) is a nonarithmetic 



probability measure on K that places no mass at — oo and has mean 



u rj(du) = E 



N 



Furthermore, fx n (dt) — rf n (dt), where rf n denotes the nth convolution of rj with itself. 

We now present a generalization of Goldie's Implicit Renewal Theorem [T3j that will 
enable the analysis of recursions on weighted branching trees. Note that except for the 
independence assumption, the random variable R and the vector (iV, C\, C2, . . . ) are 
arbitrary, and therefore the applicability of this theorem goes beyond the recursions 
that we study here. Throughout the paper we use g(x) ~ f(x) as x — > 00 to denote 
lim x - ) . 00 g(x)/f(x) = 1. 

Theorem 3.1. Let (N, C\, C2, ■ ■ ■) be a nonnegative random vector, where N E N U 
{00}. Suppose that there exists j > 1 with P(N > j,Cj > 



measure P(\og Cj € (iu, Cj > 



< E 



Eli c; log c,- 



< 00 



0,iV 



j) is nonarithmetic. 



> 



Ef=i C] 



= I, E 



< 7 < a, and that R is independent of (N, Cx^C^,- 
< P < a. If 



0) > such that the 
Assume further that 

Ej=i < 00 f or some 

. ) with E[RP] < 00 for any 



P(R >t)-E 



then 



N 



p{r > t) ~ #r 

where < H < 00 is given by 



t a ~ 1 dt < 00, 



(5) 



< — > 00, 



Eli^iogQ 



P(R >v)-E 



N 



3=1 



dv. 



Remarks: (i) As pointed out in [T3], the statement of the theorem only has content 
when R has infinite moment of order a, since otherwise the constant H is zero, (ii) 
Similarly as in |14j . this theorem can be generalized to incorporate negative weights 
{Ci} at the expense of additional technical complications. However, when the {Ci} > 
and R is real-valued, one can use exactly the same proof to derive the asymptotics of 
P(—R > t); we omit the statement here since our applications do not require it. (iii) 
When the {log Ci} are lattice valued, a similar version of the theorem can be derived by 
using the corresponding Renewal Theorem for lattice random walks, (iv) It appears, 
as noted in [14], that some of the early ideas of applying renewal theory to study the 
power tail asymptotics of autoregressive processes (perpetuities) is due to [3T] and [TS] . 
The proof given below follows the corresponding proof in [14j . 



Proof of Theorem \3.1\ Let Tc be the weighted branching tree defined by the non- 
negative vector (N, C\,Ci, . . .). For each i e A n and all k < n define \\\k = log IIi| fc ; 
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note that 11;^, is independent of N-^ k but not of Aj| s for any < s < k — 1. Also 
note that i\n = i since i € A n . Let J 7 ^, k > 1, denote the tr-algebra generated 
by {(N l ,C iiA) ,C m ,...):i€A j ,0<j<k-l}, and let T = <r(0,fi), IL_| = 1. 
Assume also that i? is independent of the entire weighted tree, Tc- Then, for any 
t e M, we can write P(i? > e') via a telescoping sum as follows (note that all the 
expectations in ([6| are finite by Markov's inequality and (11 1) 

P{R > e*) 



fc=0 



E^ 



= E* 

fe=0 



E 



53 m- 1 \ k R>e t ) 

(i\k)eA k 



E 



53 l(U ilk+1 R>e t ) 

(i|fc+l)eA fc+ i 



53 mi\nR > e') 

(i|n)eA„ 



5] l(U i]k R > e*) - E l(ni| fc C (i | fej) i? > e*) 

(i|fc)GA fc V i=l 



53 l(ni|„i2 > e<) 

(i|n)eA„ 



E * 

(i|fe)e^fc 



1(E > e*"^*) - 53 l(C (i | fcJ) i* > e 



t-Vln 



53 mi\nR > e 4 ) 

(i|n)eA„ 



Now, define the measures /x ra according to Lemma |3.1| and let 

N 



v n (dt) = Y Hk(dt), g(t) = e Qt P(R >e t )-E 



fc=0 



E 1 (^ i?> 



r(t) = e at P(R > e*) and <5„(t) = e Qt £ 



53 l(H i{n R>e t ) 

(i|n)£A„ 



(6) 



(7) 



Recall that R and (A;| fe , CVii^ i), C(i|fc,2) ■> ■ ■ • ) are independent of from where it 
follows that 



Nil 



l(R> e t - v ^)-Y, l ( C (^) R> 

3 = 1 

Then, for any f el and n £ N, 



,t-v,i 



TV 



e*(^-*) ff (t- % ) 



~w = E s 



fc=0 



53 e«^ ff (i-%) 

(i\k)£A k 



+ S„(t) = (g*v n ^)(t)+5 n (t). 
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Next, define the operator f(t) = e (* u ) f(u) du and note that 

r(t) = (g*v n -i)(t) + 5 n {t). (8) 

Now, we will show that one can let n — > oo in the preceding identity. To this end, 
let rj{du) = /ii(du), and note that by Lemma 3.1 ry(-) is a nonarithmctic probability 
measure on K that places no mass at — oo and has mean, 



/' = 



u rj{du) = E 



N 



> 0. 



Moreover, by Lemma 3.1 



k=0 



E !(%edt) 

(i\k)eA k 



= 5>* fc (*) 



(9) 



k=0 



is its renewal measure. Since /i / 0, then (|/| * v)(t) < oo for all t whenever / 
is directly Riemann integrable. By ^ we know that g € Li, so by Lemma 9.1 
from <? is directly Riemann integrable, resulting in (\g\ * v)(t) < oo for all t. 



E 



Thus, **/)(*) = S[Er=oE(i 
{g*v){t)=E 



(i\k)£A k 



g(t — < oo, which implies that 



exists and, by Fubini's theorem, 



E E 

/c=o(i|fe)eA fc 



E^ 

fc=0 



53 e a ^(*-%) 
(i|k)eA* 



lim (g * v n ){t). 



To see that S n (t) —¥ as n — > oo for all fixed t, note that from the assumptions 



< E 



Eli log Q 



< oo, £ 



.7=1 °7 



1, and i? 



< oo for some 



■*3=i J 

< 7 < a, there exists < /3 < a such that E Eli Cj 5 < 1 (by convexity). Then, 
for such /?, 



S n (t) = / e-^- u h au E 

J —OO 



53 l(n i |„i?>e 1 ' 

(i|n)£A„ 



rfu 



53 / e^ 1 1 (n 1]n iZ > e u ) 

\ A J — OO 



(i|n)GA„ 



E 

(i\n)GA n 



min{t,log(ni|„_R)} 



°/ u du 



< 



-E 



E ^\n R f 
(i|n)eA„ 



(10) 
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It remains to show that the expectation in (10 1 converges to zero as n — > oo. First 
note that from the independence of R and 7~c, 



E 



E 

(i|n)eA„ 



= E[R ]E 



E aW 

(i\n)eA n 



where E[RP] < oo, for < j3 < a. For the expectation involving ILi^ condition 
on T n -\ and use the independence of (A^i n _i, Cti\ n -i t u, G(i\ n -i^), ■ ■ • ) from T n -\ as 
follows 



E 



E 

(i|n)eA„ 



= E 



= E 



= E 



E * 

(i|n-l)e.A n _i 



E ( n ii«-i)^ 

(i|n-l)GA„_i 



AT, i 



E C (i\n-1,. 
3=1 



3) 



•7"n-l 



E^ 

3 = 1 







(■ 






3=1 



E ( n ii— i)" 3 

(i|n-l)eA„_! 



(iterating n — 1 times). (11) 



Since £7 



< 1, then the above converges to zero as n —¥ oo. Hence, the 



preceding arguments allow us to pass n — > oo in Q, and obtain 

r(t) = (g*v)(t). 



Now, by the key renewal theorem for two-sided random walks, see Theorem 4.2 in 



7j, 



' / v"P(i? > v) efo = r(i) ->• - / 5(u) du = H, t -> oo. 



Clearly, H > since the left-hand side of the preceding equation is positive, and thus, 
by Lemma 9.3 in [T4"] . 



P{R> t) ~ Ht- a , t^oo. 
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Finally, 



H = 



1 
H 
1 

n 
1 

M J- 
1 

M ■ 
1 

A* Jo 



,-(«-*) 



g(t) (it dit 



-oo ./ — oo 
oo 

t 



e*fl(t) y e- u dudt 
g{t) dt 



P(R >e t )-E 



P{R >v)-E 




4. The linear recursion: R — ^, =1 CiRi + Q 

Motivated by the information ranking problem on the Internet, e.g. Google's PageR- 
ank algorithm |191 127] , in this section we apply the implicit renewal theory for trees 
developed in the previous section to the following linear recursion: 



N 

R = CiRi 

i=l 



Q, 



(12) 



where (Q, N, C±, C2, ■ ■ ■ ) is a nonnegative random vector with N E N U {00}, 
P(Q > 0) > 0, and {i?i}i S N is a sequence of iid random variables independent of 
(Q, N, C\, C2, . . . ) having the same distribution as R. Note that the power tail of R 
in the critical homogeneous case (Q = 0) was previously studied in [53] and [T7]. In 
Section |4.3| we will give an alternative derivation of those results using our method and 
will provide pointers to the appropriate literature. 

As for the nonhomogeneous case, the first result we need to establish is the existence 



and finiteness of a solution to (12). For the purpose of existence we will provide an 



explicit construction of the solution R to ( 12 ) on a tree. Note that such constructed 



R will be the main object of study of this section. 

Recall that throughout the paper the convention is to denote the random vector 
associated to the root node by (Q, N, C\, C2, . . . ) = {Q$, iVjj, C(a t \\, C(0,2)i • • •)• 

We now define the process 



W = Q, W n = ]T QiBi, n>l, 



(13) 



on the weighted branching tree 7q,C) as constructed in Section [2] Define the process 
{i?( ra )} n >o according to 



R (n) = J2w k , n>0, 



(14) 



fc=0 
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that is, is the sum of the weights of all the nodes on the tree up to the nth 

generation. It is not hard to see that R^ satisfies the recursion 

N t N 

R (n) = E c m Ri r 1] + = E CjR { r 1] + Q> n>i, (15) 

3=1 3=1 

where {R^ 1 } are independent copies of corresponding to the tree starting 

with individual j in the first generation and ending on the nth generation; note that 
R^ = Qj- Similarly, since the tree structure repeats itself after the first generation, 
W n satisfies 

w n = E Qi^i 

N® „ 

= E c '(^ fe ) E Q(k,...,i n )Y[C(k,...,ij) 

fc=l (fc,...,i n )eA n 3=2 

JV 

E^^n-D.fe' ( 16 ) 
fe=l 



P 



where {^(n-i),^} is a sequence of iid random variables independent of (N, C\, C2, . . . ) 
and having the same distribution as W n —i- 

Next, define the random variable R according to 

00 

R^ lim iiW =y^W k , (17) 

k=Q 



where the limit is properly defined by (14) and monotonicity. Hence, it is easy to 
verify, by applying monotone convergence in (15), that R must solve 

N), N 

r = E c m) R{ r ] + = E c A oo) + Q> 

3=1 3=1 

where {-Rj^}jeN are iid, have the same distribution as R, and are independent of 
(Q,N,C 1 ,C 2 ,...). 

The derivation provided above implies in particular the existence of a solution in 
distribution to (12). Moreover, under additional technical conditions, R is the unique 
solution under iterations as we will define and show in the following section. The 
constructed R, as defined in (17), is the main object of study in the remainder of this 
section. 

4.1. Moments of W n and R 

In this section we derive estimates for the moments of W n and R. We start by 
stating a lemma about the moments of a sum of random variables. The proofs of 
Lemmas |4.1| |4.2| and |4.3| are given in Section |7.2| 
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Lemma 4.1. For any k G N U {00} let {Cj}|_ 1 be a sequence of nonnegative random 
variables and let {Yi}^ =1 be a sequence of nonnegative iid random variables, inde- 
pendent of the {Ci}, having the same distribution as Y . For (3 > 1 set p = \/3~\ G 
{2, 3, 4, . . . }, and if k — 00 assume that E*=i CiXi < 00 a - s - Then, 



E 



<(E[Y^]) 



E 



E^ 



Remark: Note that the preceding lemma does not exclude the case when 



E 



CiY 



= 00 but E 



(E-=i CiFi) -E-=i(^ 



< 00. 



We now give estimates for the /3-moments of W„ for f3 G (0, 1] and (3 > 1 in Lem- 
mas |4.2| and |4.3[ respectively. Throughout the rest of the paper define 



Pa 



E 



Ef =1 q 



for any /3 > 0, and p = p\. 
Lemma 4.2. For < /3 < 1 and a// n > 0, 

E[W£] < E\Q%-. 

Lemma 4.3. For [3 > 1 suppose E[Q P ] < 00, E fe^i ^ 
Then, there exists a constant Kp > such that for all n > 



< 00, and pV p@ < 1- 



Now we are ready to establish the finiteness of moments of the solution R given by 
(171 in Section [4] The proof of this lemma uses well known contraction arguments, 
but for completeness we provide the details below. 



Lemma 4.4. Assume that E[Q@] < 00 for some (3 > 0. In addition, suppose that 
either (i) pp < 1 ifO < /3 < 1, or (ii) (/? V pp) < 1 anc! E 



EtiOi) < 00 if p>i. 

Then, i?[i? 7 ] < 00 for all < 7 < /3, and m particular, R < 00 a.s. Moreover, if 
/3 > 1, i?*-™' —I i?, where Lp stands for convergence in (E\ ■ j^) 1 /^ norm. 



Remark: It is interesting to observe that for (3 > 1 the conditions pp < 1 and 



(Ef =1 



< 00 are consistent with Theorem 3.1 in Proposition 4 in [T7] 



and Theorem 2.1 in [23], which give the conditions for the finiteness of the /3-moment 
of the solution to the related critical (p\ = 1) homogeneous (Q = 0) equation. 



Proof. Let 



Then by Lemmas 4.2 and 4.3 



I Pp if P < 1 

\pV^, if/3>l. 



(18) 
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for some K > 0. Suppose j3 > 1, then, by monotone convergence and Minkowski's 
inequality, 



E[BP] 



< oo. 




This implies that R < oo a.s. When < (3 < 1 use the inequality (X)fe=o ^fc) — 
S/c=o f° r an y 2/« — instead of Minkowski's inequality. Furthermore, for any 

< 7 < A 



< OO. 



That 4 i? whenever /3 > 1 follows from noting that E[\R^ - R\P) 



E 



bound E[\R^ - R\P] < Kif +l /{l - rjW) 



and applying the same arguments used above to obtain the 



Next, we show that under some technical conditions, the iteration of recursion (12 1 
results in a process that converges in distribution to R for any initial condition Rq . To 
this end, consider a weighted branching tree 7q,C) as defined in Section[2] Now, define 

RZ^R^-V+WM), n>l, 
where i?( n_1 ) is given by (14), 

W n (R* ) = J2 R h u ^ (19) 
ieA n 

and {i?Q are iid copies of an initial value Rq, independent of the entire weighted 

tree 7q,c- It follows from (151 and (19) that, for n > 0, 



JY 



K+i =Y, c i R( r 1)+ Q+ w ^( R °) =e^ 



R 



(n-l) 



e ^ n W 



fc=2 



(20) 

where {-Ry™ 1 '} are independent copies of R( n_1 ) corresponding to the tree starting 
with individual j in the first generation and ending on the nth generation, and A n j is 
the set of all nodes in the (n + l)th generation that are descendants of individual j in 
the first generation. It follows that 

JV 

R n+1 = E ^3 R n,j + Q) 
3=1 



where {-R^} are the expressions inside the parenthesis in (20). Clearly, {R* n j} are iid 
copies of R* n , thus we show that R* n is equal in distribution to the process derived by 
iterating ( [l2| with an initial condition Rq . The following lemma shows that R* n =>■ R for 
any initial condition Rq satisfying a moment assumption, where denotes convergence 
in distribution. 
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Lemma 4.5. For any initial condition Rq > ; if E[Q@\, E[(Rq)P] < oo and 



E 



< 1 for some Q < f3 < 1, i/ien 



I?, 



with E[BP] < oo . Furthermore, under these assumptions, the distribution of R is the 
unique solution with finite fj-moment to recursion (12). 



Proof. Since Ry 1 ' — > R a.s., the result will follow from Slutsky's Theorem (sec 
Theorem 25.4, p. 332 in [ID] ) once we show that W„(Rq) 0. To this end, note that 
W„(i?g), as defined by (191, is the same as W n if we substitute the Qi by the i?g;. 
Then, for every e > we have that 

P(W n (R* ) >e)< e-?E[W n (R* f] 
< e^p n ,E[(R*f} 



(by Lemma 4.2 ). 



Since by assumption the right-hand side conv erges to zero as n 
R. Furthermore, E[RP] < oo by Lemma 



4.4 



oo, then R* n 
Clearly, under the assumptions, the 



distribution of R represents the unique solution to (12), since any other possible 



solution with finite /3-moment would have to converge to the same limit. 

Remarks: (i) Note that when E[N] < 1 the branching tree is a.s. finite and no 
conditions on the {Ci} are necessary for R < oo a.s. This corresponds to the second 
condition in Theorem 1 of [IT] . (ii) In view of the same theorem from [11] , one could 
possibly establish the convergence of R* n R < oo under milder conditions. However, 
since in this paper we only study the power tails of R, the assumptions of Lemma |4.5| 
are not restrictive, (iii) Note that if E Y2i=i = 1 with a £ (0, 1], then there might 



not be a < j3 < a for which E 
that was studied in [231. 



EN r ,p 
»=1 °i 



< 1, e.g., the case of deterministic Cj's 



4.2. Main result 

We now characterize the tail behavior of the distribution of the solution R to the 



nonhomogeneous equation (12), as defined by (17) 



Theorem 4.1. Let (Q, N, C\, C%, . . . ) be a nonnegative random vector, with N G N U 
{oo}, P(Q > 0) > and R be the solution to (12) given by (|17[). Suppose that there ex- 



ists j > 1 with P(N > j, Cj > 0) > such that the measure P (log Cj € du, Cj > 0, N > j) 



nonarithmetic, and that for some a > 0, E[Q a ] < oo, < E 1°§ Q 



< oo 



d E y^L] Cf — 1. In addition, 



1. E 

2. E 

Then, 



Ef =1 ci 



< 1 and E 



Ef =1 c, 



< oo, if a > 1; or, 



(j2 N 1 c a/{1+t) 



l+e 



< oo for some < e < 1, ifO < a < 1. 
P(R>t)~Ht- a , t^oo, 
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where < H < oo is given by 



H 



E 



P(R >v)-E 



N 



J2 WiR > v) 



,i=i 



dv 



E 






a -YstiiCiRiT 


aE 


EtiCfiogc; 





Remarks: (i) The nonhomogeneous equation has been previously studied for the 
special case when Q and the {Ci} are deterministic constants. In particular, Theorem 5 



of [25j analyzes the solutions to ( 12 1 when Q and the {Ci} are nonnegative deterministic 
constants, which, when J2i=i Cf = 1, a > 0, implies that C$ < 1 for all i and 
J2i Cf log C < 0, falling outside of the scope of this paper. The solutions to ( 12 ) for 
the case when Q and the Ci 's are real valued deterministic constants were analyzed in 
[5]. For the very recent work (published on arXiv after the first draft of this paper) 
that characterizes all the solutions to (12) for Q and {C} random see 0], (ii) When 

a > 1, the condition E fEi=i < oo is needed to ensure that the tail of R is 

not dominated by N. In particular, if the {Ci} are iid and independent of N, the 
condition reduces to E[N a ] < oo since -E[C Q ] < oo is implied by the other conditions; 
see Theorems 4.2 and 5.4 in [19]. Furthermore, when < a < 1 the condition 

= 1, and 



E 



(Ef=i 



< oo is redundant since E 



the additional condition E 



(Ef=i 



< s 



i=l °i 



1+e 



< oo is needed. When the {Ci} are 



iid and independent of TV, the latter condition reduces to E[N 1+e ] < oo (given the other 
assumptions), which is consistent with Theorem 4.2 in [19]. (iii) Note that the second 
expression for H is more suitable for actually computing it, especially in the case of 
a being an integer, as will be stated in the forthcoming Corollary |4.1| When a > 1 



is not an integer, we can derive an explicit upper bound on H by using Lemma 4.6 

/ k \ a k 

Regarding the lower bound, the elementary inequality [J2i=i x i) — J2i=i 
a > 1 and Xi > 0, yields 



for 



H > 



E[Q a ] 



aE 



ELcfioga 



> 0. 



Similarly, for < a < 1, using the corresponding inequality ( E»=i x * 
for < a < 1, x 4 > 0, we obtain H < E[Q a ]/(uE Yn=i C i lo gC 



< Ei=i*< Q 

(iv) Let us 

also observe that the solution R, given by |17j), to equation (12) may be a constant 
(non power law) R = r > when P(r = Q + r E^Li C) = 1. However, similarly as in 
remark (i), such a solution is excluded from the theorem since P(r = Q+r J^ili Ci) = 1 
implies £E ( Cf log CJ < 0, a > 0. 

Before proceeding with the proof of Theorem |4.1| we need the following two technical 
results; their proofs are given in Section |7.3| Lemma |4.6| below will also be used in 
subsequent sections for other recursions. With some abuse of notation, we will use 
throughout the paper maxi<i<jv Xi to denote sup 1<i<JV+1 X4 in case N = 00. 



Implicit Renewal Theory on Trees 



15 



Lemma 4.6. Suppose (tV, Ci, C%, . . . ) is a nonnegative random vector, with N € N U 
{00} and let {i?i}igN be a sequence of iid nonnegative random variables indepen- 
dent of (N, Ci, C2, ■ ■ ■ ) having the same distribution as R. For a > 0, suppose that 

< 00 for any < (3 < a. Furthermore, assume 



that E 



(t, n 1 c a/[1+c) 



< 







' N 

E 



< 00 for some < e < 1 . TTien 
1(M > t) - / ' 



max CiRi > t 

Ki<N , 



dt 



= l -E 
a 



' N 



max Cii?, 

Ki<N 



< 00. 



Lemma 4.7. Lei (Q, N, C\, C2, ■ ■ ■ ) be a nonnegative vector with N € NU{oo} otic! Zei 
be a sequence of iid random variables, independent of (Q, AT, C\, G%, . . .). Suppose 

that for some a > 1 we have E[Q a ] < 00, E I 5Z i=1 ) < i?[i?' 3 ] < 00 /o 
< /3 < a, and 53i=i Ci-Ri < 00 a - s - TTien 



/or any 



E 



E 



(CiRi) a 



< 00. 



Proof of Theorem \4-l\ By Lemma 4.4 we know that J5[iZ^] < 00 for any < (3 < a. 



To verify that E Yli=i ^7 < 00 f° r some < 7 < a note that if a > 1 we have, by 
the assumptions of the theorem and Jensen's inequality, 



E 



N 



E c ? 



< E 



N 



vi=l 



< E 



N 



7/0 



< 00 



for any 1 < 7 < a. If < a < 1, then for 7 = a(l + e/2)/(l + e) < a we have 



A? 



E^; 



«/(l+e) 



< £ 



\i=l 



a/(H 



l + e/2 
1+e 



< 00. 



The statement of the theorem with the first expression for H will follow from 
Theorem 3.1 once we prove that condition ^ holds. To this end, define 



JV 



R* — CiRi 



Then, 

P(R >t)-E 



N 



l (CiRi > t) 



< 



i=l 



P(R >t)-P\ max dR t > t 

1 Ki<N 



P max CiRi > t 

Ki<N 



E 



N 



Y J KC i R i >t) 
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Since R = R* > maxi<j<jv CiRi, the first absolute value disappears. For the second 
one, note that 



E 



N 



E 



Now it follows that 



N 

YHC i R i >t) 



- P max CiR, > t 



Ki<N 



i=l 



E 



1 I max CiRi > t 

Ki<N 



> 0. 



P(R >t)-E 



N 



J2 HCiRi > t) 



< P(R >t)-P\ max dRi > t 

\Ki<N 



E 



N 



liQRi > t) 



P 

max CiRi ^ £ 

Ki<N 



(21) 



Note that the integral corresponding to (21 ) is finite by Lemma 4.6 if we show that 
the assumptions of Lemma |4.6| are satisfied when a > 1. Note that in this case we can 
choose e > such that a/(l + e) > 1 and use the inequality 



E^ E 



(22) 



for /3 > 1, Xi > 0, k < oo to obtain 



E 



N 



y^ c p/(i+£] 



< E 



N 



< oo. 



Therefore, it only remains to show that 



J (p(R >t)-P f mw C^ >t ] j^j t a - x dt < oo. 



(23) 



v 



To see this, note that R = R* and 1(R* > t) — l(ma,xi<i< N CiR t > t) > 0, and 
thus, by Fubini's theorem, we have 



J (p(R >t)-P ^rnas C.R, > t\ J i"" 1 dt = ^E 



(R*y 



max Cii?, 

Ki<7V 



If < a < 1, we apply (22) to obtain 



(iT) a - max C,R, 

1 KKJV 



< E 



N 
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which is finite by Lemma 4.6 and the assumption E[Q a ] < oo. 

( k \ a k 

If a > 1, we have (X)i=i a; i) — i x fi x i > 0, fc < oo, implying that we can 
split the expectation as follows 



E 



(R*) a - max dR, 

1 Ki<N 



E 



N 



i=i 



+ E 



y^fjOiRi)"- ( max C.rA 



which can be done since both expressions inside the expectations on the right-hand 
side are nonnegative. The first expectation is finite by Lemma 4.7 and the second 



expectation is again finite by Lemma 4.6 



Finally, applying Theorem |3 . 1 1 gives 

P{R > t) ~ Ht-°, 

where H = (e [£f =1 Cf log C 3 ] ) ~* / °° v^ 1 (p{R > v) - E [£f =1 1(C,-JZ > v) 
To obtain the second expression for H note that 



dv. 



P{R >v)-E 



N 



dv 



v^E 



dv 



E 



-E 



a 



E Wi R > v 

\i=l I i=l 

l VQ_1 ^ ^ ^ + Q > ^) _ XJ > U ) ) dV 

)? =1 CiRi+Q N 

i=l 

N \ a N 

Y / C l R l + Q) - y2(CiRi) a 



N rCiR, 



v a ~^dv 



(24) 
(25) 



where (24 1 is justified by Fubini's Theorem and the integrability of 

/ N \ N 



< v c 



l\Y^C l R l + Q>v\ - J2 HCiRt > v) 

\i=l / i=l 

- 1 (l (j2c l R l + Q>v^J -1 ^ 
v a-i f UQRi >v)-l( max C l R l > 

X 1 — 1 \l<i<N 



max CiRi > v 

<i<N 
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which is a consequence of (23) and Lemma 4.6 and (25) follows from the observation 
that 



N 



l[Y J C i R i + Q>v 



N 

and v a - 1 ^l(C l R i > v) 
»=i 



are each almost surely absolutely integrable with respect to v as well. 
This completes the proof. 

As indicated earlier, when a > 1 is an integer, we can obtain the following explicit 
expression for H . 



Corollary 4.1. For integer a > 1, and under the same assumptions of Theorem 4-1 



the constant H can be explicitly computed as a function of E[R k ], E[C k ], E[Q k ], < 
k < a — 1 . In particular, for a = 1 . 



H 



E[Q] 



E 



and for a 



H = 



E[Q 2 ] + 2E[R]E 




c- 


+ 2(E[R]) 2 E 




2E 


E£iC?logCi' 





E[R} = 



E[Q] 



1 - B 



Ef =1 a 



Proof. The proof follows directly from multinomial expansions of the second expres- 
sion for H in Theorem 14. II 



4.3. The homogeneous recursion 

In this section we briefly describe how the methodology developed in the previous 
sections can be applied to study the critical, E Ej=i @i = 1> homogeneous linear 
recursion 



i? = CiRi, 

i=l 



(26) 



where (iV, Ci, C*2, . . . ) is a nonnegative random vector with TV € NU {oo} and {i?i}igN 
is a sequence of iid random variables independent of (N, C±, C2, ■ ■ ■ ) having the same 
distribution as R. This equation has been studied extensively in the literature under 
various different assumptions; for recent results see j23[T71l2] and the references therein. 

Based on the model from Section [4] we can construct a solution to (26) as follows. 
Consider the process {W n } n > defined by (13) with Q ; = 1. Then, the {W n } satisfy in 

distribution the homogeneous recursion in ( 16 ) and, given that E Ef=iC* = 1, 

wo 

have i?[W n ] = 1. Hence, {Wnjn^o is a nonnegative martingale and by the martingale 
convergence theorem W n — > R a.s. with E[R] < 1. Next, provided that 



E 



' N 



CilogQ 



< and E 




< 00 
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it can be shown that E[R] = 1, see Theorem 1.1(d) in [2] (see also Theorem 2 in [23"]); 
log + x — max(logx, 0). Furthermore, as argued in equation (1.9) of [2J, it can easily 
be shown that this R is a solution to ( 26 1 . Note that the same construction of the 



solution R on a branching tree was given in [2J and [23] . Since the solutions to ( 26 1 
are scale invariant, this construction also shows that for any m > there is a solution 
R with mean to; or equivalently, it is enough to study the solutions with mean 1. 
Moreover, under additional assumptions it can be shown that this constructed R is 
the only solution with mean 1, e.g. see [22J [231 [T7] . However, it is not the objective of 
this section to study the uniqueness of this solution, rather we focus on studying the 



tail behavior of any such possible solution (since our Theorem 3.1 does not require 
the uniqueness of R). As a side note, we point out that |26| ) can have solutions if 

E \Y,?=i Cf 1 = 1 for some < /3 < 1, as studied in [22JEZ]. 

A version of the following theorem, with a possibly less explicit constant, was 
previously proved in Theorem 2.2 in [23] and Proposition 7 in [17]; they also study 
the lattice case. Regarding the lattice pointed out earlier in the remark after 

Theorem |3.1| all the results in this paper can be developed for this case as well by 
using the corresponding renewal theorem. 

Theorem 4.2. Suppose that there exists j > 1 with P(N > j, Cj > 0) > such that 
the measure P(logCj € du, Cj > 0, N > j) is nonarithmetic. Suppose further that for 



some a > 1, E 



(Ef =1 



< oo, E 



Ef=iCf fog + a 



< oo and E 



Ef =1 Q 



1. Then, equation (26) has a solution R with < E[R] < oo such that 

P(R>t)~Ht- a , t^oo, 
where < H < oo is given by 
1 



H = 



E 



Eti^fiogQ 



P(R > v) — E 



N 



HCiR > v) 



dv 



E 


^EilrCA) 




aE 




iCf log C- 





Furthermore, if P(N > 2) > 0, N = J^Zi K C i > °)> then H > 



Proof. By the assumptions, the function ip(9) = E Et=i Cf 
continuous on [1, a], since <p(l) — <p(a) — 1. Furthermore, by standard arguments, 
it can be shown that both ip'(&) and f"(0) exist on the open interval (l,ct) and, in 
particular, 

N 



is convex, finite, and 



<(/' (9) = E 



Clearly, ip"(6) > provided that P(d € {0, 1}, 1 < i < N) < 1. To see that this is 



indeed the case, note that E 



C t = 1 implies that P{C l = 0, 1 < i < N) < 1, 
which combined with the nonarithmetic assumption yields P(C, € {0, 1}, 1 < i < N) < 
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1. Hence, there exists 1 < 6\ < 9 2 < a such that ip'(0i) < and ^'(#2) > 0, implying 
by the monotonicity of <p'(-) and monotone convergence that 



< <p'(a-) = E 



N 



^CflogC, 



ip'{l+) = E 



< E 



N 



N 



^Cflog+C, 



< 00 



and 



(27) 



log d 



< 0. 



The last expression and the observation E ^(j^f =1 CA log H 



Ef =1 c, 



< 00 (implied 



by E 



(EL 



< 00) yields, as argued at the beginning of this section, that 

recursion ( |26[ ) has a solution with finite positive mean, see Theorem 1.1(d) and equation 
(1.9) in gJTsee also Theorem 2 in |2"3"1) . 

Next, in order to apply Theorem 3.1 we use p7b and Els?] < 00 for any < ft < a; 
the latter follows from Theorem 3.1 in [2] and the strict convexity of ip(-) argued above 
(see also, Proposition 4 in pjj and Theorem 2.1 in [23]). The rest of the proof, except 
for the H > part, proceeds exactly as that of Theorem |4.1| by setting Q = 0. 

Regarding the H > statement, the assumption P(N > 2) > implies that there 
exist 1 < n < 00 and l<ii<i 2 <n + l such that P(N = n, C i± > 0, C i2 > 0) > 0, 
which further yields, for some S > 0, 

P{N>i 2 ,C il >S,C i2 >S) >0. (28) 

Next, by using the inequality (ari +^2)° > x" + x 2 for Xi,x 2 > and a > 1, the 
second expressions for H in the theorem can be bounded from below by 

H > E [l(N > » 2 ) {{Ci^ + C l2 R l2 ) a - (CM- ~ (G l2 R l2 r)} 

To further bound the numerator in ( |29[ ) we define the function 

f(x) = (1 + x) a - 1 - x a - cx a '\ 

where 0<e<a — 1, 0<c<2 7 — 1 and 7 = a — 1 — e. It can be shown that f(x) > 
for x € [0, 1], since /(0) = and /'(i) > ax 7 ((l + 1/x) 7 - 1 - c) > on [0, 1]. Hence, 
by using the inequality f(x) > 0, we derive for x\ > 0,x 2 > 0, m.ax.{xi,x 2 } > and 
x = min{.Ti, ar2}/max{a;i, x 2 } 

[xx + x 2 ) a -x?-x%= (max{xi, x 2 }) a ((1 + x) a - 1 - x Q ) 

> c(max{ii, X2})"a; Q_e > c(min{xi, ^2})"; 

the inequality clearly holds even if max{ii,i2} = since both of its sides are zero. 



Thus, by applying this last inequality to (|29j) and using ([28]), we obtain 
cE [1(N > i 2 ) (imn{C il R il ,C h R i2 }) a ] 



H > 



aE 



> c6 a P{N > i 2 ,C tl > S,C l2 > S)E[(mm{R tl ,R l2 }) a ] Q 



aE 



This completes the proof. 
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Remarks 
of Theorem 



by E 



4.1 



i Note that the assumptions of Theorem 4.2 differ slightly from those 
in the condition 



< E 



EtiCfiogC 

< oo, the strict convexity of ip{9) 



< 



E 



oo, which is implied 
and the 



EN n b 
i=l °J 



hypothesis that (p(l) = <p{a) — 1, as argued in the preceding proof, (ii) The assumption 
P(N > 2) > is the minimal one to ensure the existence of a nontrivial solution, see 
conditions (HO) in [22] and (C4) in [2]. Otherwise, when P(N < 1) = 1, W n is a simple 
multiplicative random walk with no branching; clearly, in this case our expression for 
H reduces to zero. Also, if P (Ei=i C< = 1) = 1, i? can only be a constant; see 
the remark above Theorem 2.1 in [23] • However, this last case is excluded from the 
theorem since P(^2 i=1 Ci — 1) — 1 implies Cj < 1 a.s., which, in conjunction with 
ip(a) = 1, a > 1, yields P(C\ € {0, 1}, 1 < i < N) = 1, but this cannot happen due 
to the nonarithmetic assumption, (iii) Note also that condition (C3) in 2 (equivalent 
to P(Ci £ {0, 1}, 1 < i < N) < 1 in our notation) is implied by the nonarithmetic 
assumption of our theorem. Interestingly enough, if this last condition fails, Lemma 1.1 
of [32] shows that equation (26) has no nontrivial solutions, (iv) As stated earlier, a 
version of this theorem was proved in Theorem 2.2 of [23], by transforming recursion 



(261 into a first order difference (autoregressive/perpetuity) equation on a different 
probability space, see Lemma 4.1 in [23] , However, it appears that the method from 
[23] does not extend to the nonhomogeneous and non-linear problems that we cover 
here, since the proof of Lemma 4.1 in [23J critically depends on having both E[R] = 1 



and E 



El, a 



= i. 



Similarly as in Corollary |4.1[ the constant H can be computed explicitly for integer 
a > 2. 



Corollary 4.2. For integer a > 2, and under the same assumptions of Theorem \4-S\ 

the constant H can be explicitly computed as a function of E[R k ], E[C k ], 1 < k < a — 1. 
In particular, for a = 2, 



H = 



E 




E 







Proof. The proof follows directly from multinomial expansions of the second expres- 
sion for H in Theorem 14.21 



We also want to point out that for non-integer general a > 1 we can use Lemma 4.1 
to obtain the following bound for H, 



H < 



(E [RP- 1 ]) a/(p - 1) E 




ot - 


aE 







where p = \a] . 
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5. The maximum recursion: R — (\J^_ 1 CiRij V Q 

In order to show the general applicability of the implicit renewal theorem, we study 
in this section the following non-linear recursion: 

R= ^\/ diuj V Q, (30) 

where (Q, N, Ci, C%, ■ ■ ■ ) is a nonnegative random vector with N E N U {oo}, 
P(Q > 0) > and {i?i}i S N is a sequence of iid random variables independent of 
(Q, N, C\, C2, . . . ) having the same distribution as R. Note that in the case of page 
ranking applications, where the {R{\ represent the ranks of the neighboring pages, 
the potential ranking algorithm defined by the preceding recursion, determines the 
rank of a page as a weighted version of the most highly ranked neighboring page. 
In other words, the highest ranked reference has the dominant impact. Similarly to 
the homogeneous linear case, this recursion was previously studied in [6] under the 
assumption that Q = 0, N — 00, and the {Ci} are real valued deterministic constants. 
The more closely related case of Q = and {Ci} > being random was studied earlier 
in |18j . Furthermore, these max- type stochastic recursions appear in a wide variety of 
applications, ranging from the average case analysis of algorithms to statistical physics; 
see [T] for a recent survey. 



Using standard arguments, we start by constructing a solution to (30) on a tree and 
then we show that this solution is finite a.s. and unique under iterations and some 
moment conditions, similar to what was done for the nonhomogeneous linear recursion 
in Section |4] Our main result of this section is stated in Theorem 15. II 

Following the same notation as in Section |4j define the process 

V n = \/ Qilli, n > 0, (31) 
ieA n 

on the weighted branching tree 7q,Co as constructed in Section [2] Recall that the 
convention is that (Q, N, C\, C2, ■ ■ ■ ) = (Q0, N$, C(0,i), C(0,2), • ■ ■ ) denotes the random 
vector corresponding to the root node. 

With a possible abuse of notation relative to Section^ define the process {R^}„>o 
according to 

n 

R {n) =\/V k , n>0. 

k=0 

Just as with the linear recursion from Section |ij it is not hard to see that R^ satisfies 
the recursion 

R (n) = ( V C i*j) R{ r l) ] V O0 = ( \/ CjRf-V ) V Q, (32) 



where {R^j 1 } are independent copies of R^ n ^ corresponding to the tree starting 
with individual j in the first generation and ending on the nth generation. One can 



Implicit Renewal Theory on Trees 



23 



also verify that 



JY 



V n — \J C(0,ife) \/ Q(k,...,i n ) JJ C{h,...,ij) — V ^fc^C™-!),*-'' 

(fc,...,i„)eA„ i=2 



fe=l 



k=l 



where {V^-i)^} is a sequence of iid random variables independent of {N, C\, C2, . . . ) 
and having the same distribution as V n -\. 

We now define the random variable R according to 

00 

R 4 lim i? (n) = W V fe . (33) 

k=0 

Note that .RW is monotone increasing sample-pathwise, so R is well defined. Also, by 
monotonicity of R^ , ( 32 ) and monotone convergence, we obtain that R solves 




v q = ( v c J R( r ] 1 v 0. 



where {i?^°°' l }j e N are iid copies of R, independent of (Q, N, C±, C2, . . . ). Clearly this 
implies that R, as defined by (33), is a solution in distribution to p0|. However, this 



solution might be 00. Now, we establish the finiteness of the moments of R, and in 
particular that R < 00 a.s., in the following lemma; its proof uses standard contraction 
arguments but is included for completeness. 

Lemma 5.1. Assume that pp = E X^=i < 1 an d E[Q^] < 00 for some j3 > 0. 
Then, i?[i? 7 ] < 00 for all < "f < (3 , and in particular, R < 00 a.s. Moreover, if 
P > 1, — I i?, where Lp stands for convergence in (E\ ■ l^) 1 /^ norm. 

Proof. By following the same steps leading to (|TT|), we obtain that for any k > 0, 



< E 



(34) 



Hence, 



k=0 



< E 



_fe=0 



< 1 < 00, 

1 - P/3 



implying that _E[i? 7 ] < 00 for all < 7 < /3. 



That flOO 4 i? whenever /3 > 1 follows from noting that E[\R^ - R\P] < 



E 



for E[V£] 



< E 



fc=n+l r fc 



and applying the preceding geometric bound 



Just as with the linear recursion from Section |4j we can define the process {R^} as 
R* n ^R^VV n (R* ), n>l, 
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where 



V n (R* ) = \/ KaKu 



(35) 



and {i?Q are iid copies of an initial value Rq, independent of the entire weighted 

tree Tq,c- It follows from (32 1 and (35) that 



N 



R {n-l) v 



V RSaUcu 



Q = V Cj'Kj V Q, 



k=2 



where {Rj } are independent copies of i?( n_1 ) corresponding to the tree starting 
with individual j in the first generation and ending on the nth generation, and A n j 
is the set of all nodes in the (n + l)th generation that are descendants of individual j 
in the first generation. Moreover, {i?* ■} are iid copies of J?*, and thus, i?* is equal 
in distribution to the process obtained by iterating ( 30 ) with an initial condition Rq . 



This process can be shown to converge in distribution to R for any initial condition Rq 
satisfying the following moment condition. 

Lemma 5.2. For any R*o > 0, if E[Q%E[{R* Y] < oo and Pf) < 1 for some f3 > 0, 
then 

R* => R, 

with E\R"\ < oo. Furthermore, under these assumptions, the distribution of R is the 



unique solution with finite [3-moment to recursion (30) 



Proof. The result will follow from Slutsky's Theorem (see Theorem 25.4, p. 332 in 
[10]) once we show that V n (R^) => 0. To this end, recall that Vn(Ro) is the same as 
V n if wc substitute the Q\ by the Rq ;. Then, for every e > we have that 

P(V n (R* Q ) >e)< e-^E[V n {R* f] 

< e^p n p E[(R* f] (by @). 



Since by assumption the right-hand side conv erges to zero as n 
R. Furthermore, E[R@] < oo by Lemma 



5.1 



oo, then R* n ^ 
Clearly, under the assumptions, the 



distribution of R represents the unique solution to (301, since any other possible 



solution with finite /3-moment would have to converge to the same limit. 

Now we are ready to state the main result of this section. 

Theorem 5.1. Let (Q, N, C\, C2, . . . ) be a nonnegative random vector, with ]V£NU 
{00}, P(Q > 0) > and R be the solution to (30) given by (33). Suppose that 



there exists j > 1 with P(N > j, Cj > 0) > such that the measure P(log Cj G 
du,Cj > 0, N > j) is nonarithmetic, and that for some a > 0, E[Q a ] < 00, < 



E ^2n—xCf\ogCi < 00 and E EiLi Cf —\. In addition, assume 



N 



1. E 



2. E 



(Eilra 

f J2 N 1 c a/{1+c) 



< oo, , if a > 1; or, 

l+£" 



< 00 for some 0<e<l, i/0 < a < 1. 
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Then, 

P(R>t)~Ht- a 1 t^oo, 
where < H < oo is given by 



H = 



1 



E 



Eti^fiogQ 



P(R >v)-E 



N 



dv 



E 






Q vQ Q -Ef=i(W 


aE 







Proof. By Lemma 



■5.1 



we 



know that E\RP] < oo for any < /3 < a. The same 



arguments used in the proof of Theorem 



4.1 



give that E 



< oo for some 



< 7 < a. The statement of the theorem with the first expression for H will follow 
from Theorem 3.1 once we prove that condition ([5| holds. Define 



R* = ( V CiRij V Q. 



Then, 

P(R >t)-E 



N 



J2 WiRi > t) 



.i=i 



< 



P(R >t)-P[ max CiR z > t 

1 Ki<N 



p max CiRi > t 

Ki<N 



E 



N 



KCiRi > t) 



Since R = R* > mssL\<i<N CiRi, the first absolute value disappears. The integral 
corresponding to the second term is finite by Lemma |4.6| just as in the proof of 
Theorem 4.1 To see that the integral corresponding to the first term, 

J (p(R >t)-P ^max^ dRi > tYj t^ 1 dt, 

is finite we proceed as in the proof of Theorem 4.1 First we use Fubini's Theorem to 
obtain that 

J (p(R >t)-P ^max^ C i R l > tj \ t a ~ x dt 



-E 



< 



a 



l -E 
a 

E[Q C 
a 



(R* 



max CiR, 

Ki<N 



N 



N 



\JCiRi) V Q a - \/ dRi 



\i=l 
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Now, applying Theorem |3 . 1 1 gives 

P(R > t) ~ Ht- a , 



where H = (e [^f =1 Cf log Cj] ) * / °° u"" 1 (p(R > v) - E [^f =1 > «) 

The same steps used in the proof of Theorem |4.1| give the second expression for H 



dv. 



6. Other recursions and concluding remarks 

As an illustration of the generality of the methodology presented in this paper, we 
mention in this section other recursions that fall within its scope. One example that 
is closely related to the recursions from Sections [4] and [5] is the following 



N 



R = \/ dRi + Q, 



(36) 



\i=l 



where (Q, N, Ci, C 2 , . . . ) is a nonnegative vector with N e N U {oo}, P(Q > 0) > 0, 
and is a sequence of iid random variables independent of (Q, N, C\, C%, . . .) 

having the same distribution as R. Recursion ( 36 1 was termed "discounted tree sums" 



in pQ ; for additional details on the existence and uniqueness of its solution see Section 
4.4 in [J. 

Similarly as in [14,, it appears that one could study other non-linear recursions on 
trees using implicit renewal theory. For example, one could analyze the solution to the 
equation 



CiRi + E>i\J R, 



Q, 



where (Q, N, d, C 2 , . . . ) is a nonnegative vector with N e N U {oo}, P(Q > 0) > 0, 
and {i?, Ri}i>i is a sequence of iid random variables independent of (Q, N, Ci, C%, . . . ). 
Here, the primary difficulty would be in establishing the existence and uniqueness of 
the solution as well as the finitcncss of moments. 



7. Proofs 

7.1. Implicit renewal theorem on trees 

We give in this section the proof of Lemma |3.1 



Proof of Lemma \3.1\ Observe that the measure E YliLi l(l°gC< G du, Ci > 0) is 
nonarithmetic (nonlattice) by our assumption since, if we assume to the contrary that 
it is lattice on a lattice set L, then on the complement L c of this set we have 



= E 



N 



^i(io g c i ei c ,c i >o) 



> P(log Cj e L c , Cj > 0, N > j) > 0, 



which is a contradiction. Hence, r\ is nonarithmetic as well, and it places no mass at 
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-oo due to the exponential term e au . To see that it is a probability measure note that 



r](du) 



e au E 



N 



E 



E 



^lOogQedu) 

3=1 

/OO 
e a "l(logC,- £du) 

N 



A'* 



(by Fubini's Theorem) 



3=1 



= 1. 



Similarly, its mean is given by 



urj{du) = E 



A 



3=1 



To show that fi„ — rj* n we proceed by induction. Let T n denote the cr-algebra 
generated by {(Ni, Cni^Cnp), ■ ■ ■ ) : i £ Ay, Q < j < n — l}, .Fo = cr(0, O), and for 
each i £ A n set V\ = loglli. Hence, using this notation we derive 



[X n+ i((-oc,t]) 



e au E 



e au E 



5^ 5^i(H+iogc (iJ) edu) 



Ai 



^l^i+logC^) £ du) 



3=1 



^l(logC (iij) Gd«-V5) 
3=1 



J' n 



E 



oo 

?/((-oo,t- w])^„(dv), 



where in the fourth equality we used the independence of (iVi, C(i,i) , C(i.2), ■ ■ ■ ) from 
Therefore /i„+i(di) = (77 * fi n )(dt) and the induction hypothesis gives the result. 

7.2. Moments of W n 

In this section we prove Lemmas |4.1[ |4.2| and |4.3| We also include a result that 



provides bounds for for integer p, which will be used in the proof of Lemma 4.3 



Proof of Lemma gj} Let p = \0] £ {2,3,...} and 7 = /3/p £ (fi/(0 + 1),1]. 
Suppose first that k £ N and define Ap(fc) = {(ji,...,jk) £ N fe : ji + ■ ■ ■ + jk — 
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P, < ji < p}- Then, for any sequence of nonnegative numbers {yi}i>i we have 



E^f+ E 

* =1 (Ju---,jk)eA p (k) 



i=l 



P 

ji, ■ • • 



yf • • • J/J* 



(37) 



E 

v (ji,---jfc)eAp(fe) 

where for the last step we used the well known inequality fj3i=i ^i) — Ei=i x l f° r 
< 7 < 1 and > 0. We now use the conditional Jensen's inequality to obtain 

f> k 



E 



E^O -E 



si=l 



i=l 



< E 



< E 



E 

v(ii,-,jfc)e^ y (fc) 



Ji, ■ ■ ■ ,3k 



(dY^ ■ ■ ■ (C k Y k y* 



(by (371) 



E 



E 

Ui,...,j k )eA p (k) 



. ' . (CiFi)^ • • • (C k Y k y* 
Ji, ■ ■ ■ , 3k 



Cx,...,C k 



E 



P 

3i, ---,3k 



, )Ct---CtE 



y3l _ _ _ yjfc 



Ci , . . . , C k 



E 

K Ui,----Jk)eA p (k) 

Since {Yi} is a sequence of iid random variables having the same distribution as Y, 
independent of the CVs we have that 



E 



C\, . . . ,c k 



E 



where ||Y|| K = (£[|Y | k ]) 1/k for k > 1 and ||Y|| = 1. Since \\Y\\ K is increasing for 
k > 1 it follows that ||F||| < 11^11^!. Hence 



\\Y\\f l ---\\Y\\f k <\\Y\\ p p _ 1 , 



which in turn implies that 



E 



E CiYi ) -£(c^ 



< E 



E 

l(ilr-j»)^((*) 



<ll^ll^ 



P 

3i, ■ • ■ , jfc 

P fc 



^■••^imij-x 



E c 



Implicit Renewal Theory on Trees 



29 



This completes the proof for k finite. 

When k — oo, first note that from the well known inequality (x% + x 2 )^ > x± + x 2 
for X\, X2 > and /3 > 1 we obtain the monotonicity in fc of the following difference 



fk+i 







fc+i 







Hence, 



E 



lim 

k— >oo 



< lim # 

k— >oo 



E^ -E( c ^: 



E 



< lim Hrii^^ 



i=l 

ii , ■ ■ ■ , h 





(38) 



(c^yi ■ ■ ■ (c k Y k y« 



\Y\\Z-iE 



E^ 



(39) 



where (38) and (39) are justified by monotone convergence. 



Proof of Lemma 4-2 



< p < 1, yi > 0, k < oo, to obtain 



We use the well known inequality Vij — Ei=i Hi f° r 



E\WS\ = E 



< E 



N 



(n-l),i 



N 



EGf<-r 



= E\W^_-^pp (by conditioning on iV, Cj and Fubini's theorem) 
< p^E[Wq] (iterating n times) 

= p n p E[Q% (40) 

Before proving the moment inequality for general /3 > 1, we will show first the 
corresponding result for integer moments. 



Lemma 7.1. Let p G {2, 3, . . . } and recall that p p = E J^iLi Cf 



, p = pi. Suppose 



E[Q P ] < oo, E (EiliCi) < o°. and pV Pp < 1- T ^en ; i/i ere exists a constant 
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K p > such that 

E[W%]<K p (pV Pp ) n 

for all n > 0. 

Proof. We will give an induction proof in p. For p = 2we have 
E\Wl\ = E 



f N \ • 



E 



N 



i=l 



A' 



E^ 



<p 2 £[T^_ 1 ]+ J E 



at 



(£;[w„_!]) 2 £; E Ci ^ 

(by conditioning on N. Ci and Fubini's theorem) 

21 



E^ 



{E[W n ^]f 



Using the preceding recursion and noting that, 

E[W n - l ] = p n - 1 E[Ql 

we obtain 



E[W 2 ] < P2 E[Wl_ 1 ]+Kp 2 ^- 1 \ 



(41) 



where K = E 



(Ef =1 



a 



(E[Q}) . Now, iterating (41) gives 
E[W 2 n ] < p 2 (p 2 E[W 2 _ 2 ] + Kp 2 ^) + Kp 2 ^ 

n— 2 

< pT 1 (p 2 E[W 2 } +K)+KJ2 A P 2(n '^ l) 

i=0 

n-1 

= ^ J B[Q 2 ] + ir^p|p 2 ("- 1 - i ) 

< (p 2 V p)"£[Q 2 ] + K( P2 V p)" E(^ V PT~ 
K - 



4=0 



< | E[Q 2 ] + V P y (p 2 V p)* 

p 2 V p ^ » 



= X 2 (P2 Vp) f 

which completes the case p — 2. 
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Suppose now that there exists a constant K„_i > such that 

Em' 1 } < (Pp-i V pf 



(42) 



for all n > 0. Then, by Lemmas 4.1 and 4.2 we have 



= E 



N 



N 



< 



p/(p-i) 



i=l 
'TV \ Pi 



AT 



V P/(P-1) 



/ N 

E^ 

\i=l 



E^-D, 

P P £ K-l] 



< £ 



N \ P 



E^ 

\i=l 



(Ap _ l) p/(p-D (pp _ 1 V p) (n-l)p/(p-D + ftJ^W^], 



where in the second equality we conditioned on N, Ci and used Fubini's theorem, 
and the last inequality corresponds to the induction hypothesis. We then obtain the 
recursion 

« < PpE\K-i] + K (Pp-i VP)^, (43) 
(Eili Ci)* (-Kp-i) p/(p_1) - Iterating @ as for the case p = 2 gives 



where K = E 



E[W%] < p n p E[QP] +kJ2p^ (Pp-i V p)^" 1 -^- 1 ) 

i=0 

71-1 

< (ft, Vp)"£[Q p ] + A'^(p p V p)« B - 1 )H/(P-i) (44) 

i=0 

n-1 

= (ft V pTeiq?} + K(p p V p)- 1 £(p p V p)("-*-i)/(P-D 

i=0 

< [e\Q p ] + K( Pp V p)- 1 f> p V p)^r J ( Pp V p)" 
= Ap(p p V p) n , 



where in (44) we used the convexity of (p{/3) — pp, i.e. p p _i = <p(p — 1) < <p(l) Vy(p) 
pV pp. 



The proof for the general /3-moment, /3 > 1, is given below. 
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Proof of Lemma \4-3\ Set p = \0\ > j3 > 1. Then, by Lemmas 4.1 and 4.2 
E[Wg] = E 



N 



\P N 



vi=l 



,/3/(p-1) 



i=l 



+ E 



N 



E 



\i=l 
t N 

V.i=l 



1=1 

E^f<-D, 



By Lemma [7TT 



N 



= p E[W^_ 1 ]+K(p p _ 1 \/p)^- 1 ^, 
where J = /3/(p — 1) > 1. Finally, iterating the preceding bound n — 1 times gives 



£7[W£] < ] + K £ Pp(P V Pp-i)^"- 1 "^ 

i=0 

< ](p V p P Y + X ]T(p V p p )^ n - l -V +i (by convexity of tp(fi) = p p ) 



i=0 



n-l 



= E[QP]{p V ppT + K(p V ppY- 1 £(p V PP) {1 ~ 1)1 
<Kp(pVp p ) n . 
This completes the proof. 

7.3. Linear nonhomogeneous recursion 

In this section we give the proofs of the technical Lemmas |4.6| and |4.7| for the linear 
recursion. 



Proof of Lemma 4-6 Note that the integral is positive since 

" N 

1 {CiRi > t) 



P [ max CilL >t \ =E 

Ki<N 



1 ( max CiRi > t 

Ki<N 



< E 



To see that the integral is equal to the expectation involving the a-moments note 
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that 



E 



E 



N 



l(CiRi > t) 



=1 

' N 



- P { max dR t > t ) t^ 1 dt 



Ki<N 



t"- 1 dt 



= E 



E 



^ l(CiRi > t) - 1 I max QRi > t 

i=l \l<i<N 

J {Y1 l ( C i R i > t )- 1 (^ 1 m .ax r CiRi >t)\ t"- 1 dt 

" N / 

V -(C l R l ) a --( max C l R l 
^ a a \i<i<N 



(by Fubini's theorem) 



where the last equality is justified by the assumption that X^i(Ci^)" < 00 a - s - 

It now remains to show that the integral (expectation) is finite. To do this let 
X = (N, C\ , Ci , . . . ). Similar arguments to those used above give 



E 



N 



JO 



= E 



E 



.4=1 



E 



P[ max C\R l >t) \t a - 1 dt 

Ki<N 



> t) 

N 

V l(C l R l >t)-li max C l R l > t 

f-^ \l<i<N 



i=l 



where in the last step we used Fubini's theorem. Furthermore, 



X 



X 



e- 1 dt 



t a - x dt 



E 



-»=i 

= £ 



V l(C 4 i?, > t) - 1 ( max Qi^ > t* 

izi v<*< N j 

N 

-1 + ^£[1(CA >t)\X] 



1 ( max Cjiij < t 

Ki<JV 



Note that given X, the random variables CiRi are independent (since the i?'s are), 
so if we let F(t) = P(R>t), then 

N N 



E 



1 f max CiRi < i 

Ki<N 



l[E[l(C i R i <t)\X} = l[(l- F(t/d)) . 

i=l i=l 



We now use the inequality 1 — x < e x for x > to obtain 

N 

Y[(l-F(t/Ci)) <e-X£int/c<). 



Next, let S — ae/ (1 + e) and set f3 = a — S. By Markov's inequality, 

N N N 

i—1 i=l i—1 
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Now, define the function g(x) = e x — 1 + x and note that g(x) is increasing for 
x > 0. Therefore, 

g(j£Fm)\ <g(t-'E[RP]£cA. 
This observation combined with the previous derivations gives 

N 



poo / 

IE l ( C i R i > t) - 1 ( max CiRi > t 

J0 [_i=l \l<t<N 



X 



t a - Y dt 



where Sp — J2^ =1 C^ and r = E[R^\ < oo. Hence, using the change of variables 
u — rSpt~@ gives 



poo 

1 + rS^) t a ~ l dt = 0r l (r8p) a / p I (e~ u -l + u) u^'^ 1 du 



Our choice of /9 = a — 6 guarantees that l<a//3 = l + e<2. To see that the 
(non-random) integral is finite note that e~ x — 1 + x < x 2 /2 and e~ x — 1 < for any 
x > 0, implying 



1 f 1 

(e" u -l + u) u- alf} - x du<- u x - alfi du 



1 



1 



2(2 - a//3) a/ (3 - 1 



*^du 



= K(3 < oo. 



Now, it follows that 



E 



N 



^2 WiRi > t) 



i=l 



P max QRi > t) ) t^ 1 dt 



Ki<N 



< E 



K(rSp) a/p ] = Kr a ' p E 



The last expectation is finite by assumption (a//3 = 1 + e), which completes the proof. 



Proof of Lemma \4- r l\ Let S — J^j— x CjRj < oo a.s., p — \a] and note that 1 < 
p-Ka. Then, since {S + Q) a - S a > and S a - ^f =1 (C,]J,) a > 0, we can break 
the expectation as follows 



E 



N 



{s + QT 



i=l 



B [(5 + Q) a — S a ] + E 



N 



N 



J2 C ^) -E( c ^) a 



i=l 



<E[(S + Q) a -S a ] + (E[R P 



-Il\ a /(P-!) 



N 



V.i=l 
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where the inequality is justified by Lemma |4.1| The second expectation is finite since 
by assumption E[RP\ < oo for any < f3 < a. For the first expectation we use the 
inequality 



(x + t) K < 



x K +t K , 



< K < 1, 



X K + K(X + t) K ~ 1 t, K>1, 



for any x,t > 0. We apply the second inequality p — 1 times and then the first one to 
obtain 

p-2 

(x + t) a <x a + a(x + t) a -H < ■■■ <x a aV* - **' + a p ~\x + t ) a - p +H p - 1 

i=i 

p-i 



<x a + a p t a +a p J2 xa ~ if 



Hence, it follows that 



p-i 



E [(5 + Q) a - S a ] < a p E[Q a ] + a p J2 E[S a - l Q 1 ]. 



(45) 



To see that each of the expectations involving a product of S and Q is finite let 
X = (Q, N, C\ , C 2 , . . . ) and note that for i = p — 1, 



E[S a - p+1 Q p ~ 1 } 



(46) 



= E 



< E 



E 



N 



a— p+1 



Q (p-l)/( a - p +l)^ c . Rj 
3=1 



iV 



)(p-l)/(o- 



- p+1) E^ 



= (£[i?])"- p+1 £: 



3=1 



N 



X 



(by Jensen's inequality) 



a— p+1 



3=1 



(47) 



where the last equality was obtained by using the independence of {Rj} and X. 

For 1 < i < p — 2 let Qi = \a — i~\ and condition on Q and X, respectively, to obtain 



E^-'Q 1 ] = E 



E 



N 



s a - 1 -J2( c 3 R ]) a ~ l Q l 



3=1 



E 



N 



3=1 



JV 



Q l E 



sr*-* - y\{c j R j ) a - i 



Q 



E[R a - l ]E 
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< E 



Q l (E[R qi ~ x \($\) E 



E[R a -' l }E 




(48) 




E[R q > 



9i-ll 



- 1 +E[R°- i } E 




(49) 



where for the inequality we used Lemma 4.1 (a — i > 1) and the inequality 5Z i= i Vi — 

/ k \P 

I X)i=i !/i ) f° r an y /3 > 1 and j/j > 0. Note that all the expectations involving R in (47 1 
and (49) are finite since E[R@] < oo for all < (i < a by assumption. Next, observe 

\ ■ f n \ "^^1 

that all the other expectations are of the form E Q l \ . =1 Cj J for 1 < i < p— 1. 

To see that these are finite use Holder's inequality with q = a /(a — i) and r = a/i to 
obtain 








< 






\i= 










= 


(■ 











IIQ'llr 



1/9 



< OO. 
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